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Studying the precision of ray tracing techniques with Szekeres models 
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The simplest standard ray tracing scheme employing the Born and Limber approximations and 
neglecting lens-lens coupling is used for computing the convergence along individual rays in mock 
N-body data based on Szekeres swiss cheese and onion models. The results are compared with the 
exact convergence computed using the exact Szekeres metric combined with the Sachs formalism. 
A comparison is also made with an extension of the simple ray tracing scheme which includes 
the Doppler convergence. The exact convergence is reproduced very precisely as the sum of the 
gravitational and Doppler convergences along rays in Lemaitre-Tolman-Bondi models. This is not 
the case when the models are based on non-symmetric Szekeres models. For such models, there is 
a signihcant deviation between the exact and ray traced paths and hence also the corresponding 
convergences. ^ 


PACS numbers: 98.80.-k, 98.80.Jk, 98.80.Es 

I. INTRODUCTION 

With the possible near-future exception of observa¬ 
tions based on gravitational waves, all astrophysical 
observations are based on light. This makes under¬ 
standing light propagation a crucial element of both 
theoretical and observational cosmology. The light that 
is observed in astrophysical observations has propagated 
through vast regions of the Universe to reach us. During 
its propagation, the light feels the exact, local spacetime 
and not some ’’average” spacetime described by a 
Friedmann-Lemaitre-Robertson-Walker (FLRW) model. 
What effects the local inhomogeneities of spacetime have 
on light propagation and how important these effects are 
is still up for debate. For redshift-distance relations it 
has however been shown that averaging over many light 
rays will reduce observations approximately to what one 
would see if the light had simply traveled through the 
averaged universe model (see e.g. [1-9]). For models 
with vanishing backreaction this implies that averaging 
over many light rays will yield results approximately 
corresponding to FLRW results. (For effects of non¬ 
vanishing backreaction on light propagation, see e.g. 
[9]. See e.g. [10-13] for introductions and reviews on 
cosmic backreaction.) The number of geodesics needed 
to obtain such results can be quite large though. In 
addition, some observables, such as reduced shear and 
CMB temperature fluctuations, are not suitable for the 
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needed averaging. It is thus important to study the 
effects of inhomogeneities on light propagation so that 
the gained knowledge can be used when interpreting 
especially high precision observations. 

Much work has gone into using perturbation the¬ 
ory to study the effects of inhomogeneities on leasing 
and redshift-distance relations (see e.g. [14-20] for 
some recent examples). Another approach is to use 
exact, inhomogeneous solutions to Einstein’s equations. 
Among the most realistic, exact solutions to Einstein’s 
equations which contain dynamical structures are the 
quasi-spherical Szekeres models [21] including their 
spherically symmetric limit, the Lemaitre-Tolman-Bondi 
(LTB) models [22-24]. Light propagation through these 
models has been vastly studied, especially with the 
purpose of studying the effects of inhomogeneities on 
CMB and supernova observations (see e.g. [7-9, 25-59] 
for some examples) With the exception of onion 
models (see e.g. [38]), these models are double or 
triple structure models (see e.g. [62]) and are thus 
not individually useful for realistic studies of light 
propagation over large distances. A possible method for 
overcoming this issue is to combine few-structure models 
to build multiple-structure swiss cheese models (first 
introduced in [63]). These models have a high degree 
of complexity and are thus very important and useful. 
However, they suffer from simplicities such as a lack of 
interaction between individual structures (’’holes”) and 
often the holes in the cheese are made up of only a few 
specific inhomogeneous models representing structure 
formation on a limited scale interval. Perhaps because 
Swiss cheese models have these insufficiencies, universe 


^ Other models breaking either the assumption of homogeneity or 
isotropy have also been used to study light propagation. Ex¬ 
amples are [60, 61] concerning light propagation in Bianchi and 
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the Szekeres models, as the Szekeres structures are more compa¬ 
rable to structures in typical N-body simulations. 
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models based on the output from Newtonian N-body 
simulations are by many considered the most realistic 
models of the real universe - despite their lack of 
relativistic effects such as cosmic backreaction [64, 65]. 
The standard cosmological setting for studying light 
propagation through an inhomogeneous universe is thus 
the output from Newtonian N-body simulations. The 
caveat is that Newtonian theory does not include a 
metric. This is problematic as the metric is crucial 
for studying relativistic light propagation and weak 
lensing. To circumvent the problem, it is standard to 
introduce the perturbed FLRW metric in the Newtonian 
gauge and use this metric to obtain approximation 
schemes for tracing rays through N-body simulations 
and computing the corresponding shear, convergence 
etc.. Work {e.g. [66-68]) has been done with the goal 
to improve the standard ray tracing schemes which 
are based on several approximations besides the use of 
the perturbed FLRW metric. However, it is difficult 
to obtain exact quantifications of the impact of these 
approximations and thus to determine to what extent 
the suggested improvements are actually worth their 
trouble. 

Exact, inhomogeneous solutions to Einstein’s held 
equations, such as the Szekeres models, are excellent 
tools for studying the validity and precision of results 
obtained by using mainstream approximation schemes. 
Such models have for instance been used to study 
the possibility of a giant void being the reason for 
the CMB cold spot [69-71], to study the precision of 
single structure thin lensing [55], and to illustrate the 
dominance of Doppler convergence at low redshifts [72] . 

In this work, quasi-spherical Szekeres swiss cheese 
and onion models will be used to study the precision of 
the simplest, standard ray tracing scheme used for ob¬ 
taining the convergence along individual geodesics. The 
convergence is the most basic cosmologically interesting 
quantity obtainable from ray tracing and can be used 
to compute e.g. distance measures. The gravitational 
convergence can also be used to obtain the shear (see 
e.g. [73]), so once the convergence has been obtained, 
many observables can be studied (see e.g. [74, 75] for 
some examples). Hence, obtaining the convergence along 
light rays is the typical goal of ray tracing. 

The precision of the convergence obtained from 
standard ray tracing schemes will here be studied by 
using the direct comparison approach introduced in 
[76]. Using this approach, Szekeres models are repro¬ 
duced as (mock) N-body data. Standard ray tracing 
approximation schemes can then be combined with 
the N-body data to obtain the convergence along light 
rays. The approximate results hereby obtained are 
compared to the exact results obtained by studying 
exact light propagation in the exact Szekeres spacetime. 
The comparisons made in the following show that in 
the special case of LTB single void/swiss cheese models, 
the ray tracing method reproduces exact ray paths and 


corresponding convergences very well. However, in the 
more general case where the single void models are 
non-symmetric, as well as in the case of onion models, 
the ray tracing scheme does not reproduce the exact 
results. 


II. QUASI-SPHERICAL SZEKERES MODELS 

This section gives a brief introduction to the quasi- 
spherical Szekeres models including a description of the 
particular models used in this work. See e.g. [62, 77] for 
more details on the Szekeres models and the procedure 
for constructing particular models. 


The quasi-spherical Szekeres models are a family of 
exact, inhomogeneous solutions to Einstein’s equations. 
Using stereographic coordinates, the spacetime of the 
model can be described by the line element: 


ds^ = —c^dV + 


$,(t,r)-$(t,r)§^|^) 


1 — k{r) 
$(t,r)^ 
E{r,p, qf 


■dr^ 


( 1 ) 


{dp^ +dq^) 


Subscripted commas followed by one or more coordinates 
or indices indicate partial derivatives with respect to the 
implied coordinate(s). 

Inserting the line element into Einstein’s equations for 
a dust universe with a cosmological constant yields the 
following two useful relations: 



2^ , 1 A ^2 

4) 3c2 


( 2 ) 


2M^r - 


/3 = SttGat/c^ 


( 3 ) 


The function M = M{r) appearing in these equations is 
a temporal integration constant depending on the radial 
coordinate and corresponds to the effective gravitational 
mass at comoving radial coordinate r. 

E is given by E = ^(p^ + 9^) - f 
where S', P and Q are continuous but otherwise arbitrary 
functions of r. The quasi-spherical Szekeres models 
reduce to LTB models when P, Q and S are constant 
functions. 


The LTB model has two free functions and an ex¬ 
tra degree of freedom from the covariance of the radial 
coordinate. A quasi-spherical Szekeres model has the 
three dipole functions as extra free functions and are 
specified by choosing dipole functions to remove the 
spherical symmetry of an underlying LTB model. 

The models studied here all have a vanishing cos¬ 
mological constant. The coordinate covariance of r is 
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FIG. 1. Present time ID density profiles of LTB and non-symmetric Szekeres models with rb = 60Mpc and varied m. For the 
non-symmetric models, the profiles are shown along the direction of highest degree of anisotropy. The comoving coordinates 
are normalized at present time in units of O.lMpc. 


in the models removed by setting = r, where 

to is the present time according to an Einstein-de 
Sitter (EdS) model with reduced Hubble parameter 
h := iJo/(100km/s/Mpc) = 0.7. One of the free func¬ 
tions of the LTB models is chosen to be the time of the 
Big Bang which is set equal to zero, i.e. all the studied 
models have $(t = 0,r) = 0. The last free function of 
the LTB models is chosen to be k(r) which is specihed 
below for the different models which will be studied. 


A. Single void/Swiss cheese models 

The Swiss cheese models studied here are based on ei¬ 
ther LTB or non-symmetric Szekeres single void models. 
The final specification of the LTB single void models is 
made by defining k{r) as: 

-r"fcmax((^) -l) ifr<rfc (4) 

0 otherwise 

The models studied here have fcmax = 8 • 10“^°. The 
parameters m and ri, are varied in order to study the 
importance of void size and shape. The specific values of 
m and ri, used here are 2,4, 6 and 30Mpc, 60Mpc, lOOMpc 
respectively. Present time density profiles of the models 
with m = 60Mpc are shown in the top row of figure 1. 
To construct non-symmetric Szekeres models, the LTB 


models are modified by dehning ^ Q = — 1) with 

n = 2.3 and a = ID present time density prohles 

of the non-symmetric models are shown in the bottom 
row of figure 1. 

B. Onion models 

Aside from swiss cheese models, also onion models 
(see e.g. [38]) will be used. Onion models have an 
oscillating density field in the radial direction. In the 
models studied here, the oscillations are damped and 
will be referred to as truncated onion models. The 
damping is achieved by combining k{r) with an expo¬ 
nentially decaying function as shown in equation (5). 
The radial coordinate where the density fluctuations 
become negligible will be referred to as the truncation 
radius. 

The final specification of the LTB onion models 
studied here is made by specifying fc(r) as: 

k{r) = -femax^^cos^ e““^ (5) 


^ The constant term, —n, in Q has no physical significance but 
is included in order to simplify a technical detail regarding ray 
tracing through the origin of single void models. 
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Density profile of onion inoclcl Olitb 


Density profile of onion model 02^!, 


Density profile of onion model OS/it 



(a)model Olitb 


(b)model 02itb 


(c)model 03itb 



comoving y-coordinate (x = 0 = 2 ) 


comoving y-coordinate (x = 0 = 2 ) 


comoving y-coordinate (x = 0 = 2 ) 


(d)model Olsz 


(e)model 02be 


(f)model 03az 


FIG. 2. Present time ID density profiles of the onion models. The comoving coordinates are normalized at present time in 
units of O.lMpc. The non-symmetric models’ density profiles are shown along the axis of largest anisotropy. 


Model 

OL 

n G 

^Isz/ltb 

0.5 

1 0.0002 

02s2/ltb 

1 

2 0.0002 

^3s2/ltb 

2 

2.3 0.001 


TABLE I. Specification of non-symmetric Szekeres onion 
models and their underlying LTB models. The model names 
are subscripted with ”sz” or ”ltb” in order to distinguish be¬ 
tween the non-symmetric and spherically symmetric versions 
of the models, a alone specifies the LTB models, while n and 
a refer to the dipole function Q used to remove their spherical 
symmetry. 

All models will have /cmax = 6.2 • 10“^° and rj, = 200Mpc. 
The parameter a is varied. For a — 0, the model is a 
true onion model with an oscillating density field for all 
values of the radial coordinate. In this particular true 
onion model, the density becomes negative and hence 
that model will not be used. Truncating the model with 
a > 0.25 removes the occurrence of negative densities. 
The density fields of the onion models used here are 
shown in figure 2 and specified in table I. The spheri¬ 
cal symmetry is removed by using the dipole function 
Q = with the values of n and cr specified in table 

I. 


density fields that look very realistic since the inhomo¬ 
geneities are arranged in voids and surrounding clusters 
whereas the LTB onion models at best can be consid¬ 
ered as models with a density distribution that has been 
averaged in angular directions - a feature which is only 
barely removed by adding anisotropy through the dipole 
functions. On the other hand, the individual single void 
structures of a swiss cheese model are restricted from in¬ 
teracting with each other, while the inhomogeneities of 
onion models have coupled dynamics. 

Single void and swiss cheese models have intrinsic back¬ 
grounds, namely the backgrounds they reduce to outside 
their structures. Onion models do not have this feature. 
This is inconvenient since an FLRW background must 
be defined in order to use standard ray tracing schemes. 
It is also an advantage of the onion models since the 
real universe does not have a background either. Newto¬ 
nian N-body simulations have implicit backgrounds that 
determine the overall dynamics of the simulation box. 
It is not clear to what extent these backgrounds affect 
the structure formation of N-body simulations. How¬ 
ever, structures in typical N-body simulations are not 
separated by extended patches where the density corre¬ 
sponds to the background value. Thus, the structures in 
N-body simulations interact freely with each other simi¬ 
lar to in onion models and unlike in swiss-cheese models. 


C. Swiss cheese vs. onion models 

As semi-realistic models of large scale structure forma¬ 
tion, the Swiss cheese and onion models each have their 
advantages and disadvantages. Swiss cheese models have 


III. CONSTRUCTING MOCK N-BODY DATA 

In order to compare ray tracing through N-body 
data with exact light propagation, the Szekeres/LTB 
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models are reproduced as mock N-body data. A detailed 
description of constructing mock N-body data from 
Szekeres models was given in [76] (see also [78]) and only 
a brief review will be given below. 

In this section, subscripts will be used to distinguish 
between coordinates in the exact Szekeres spacetime 
and the underlying FLRW spacetime of the mock 
N-body data. In all other sections, these subscripts are 
suppressed as it should be clear from the context which 
spacetime the coordinates refer to. It is emphasized 
that the exact models and their N-body counterparts 
do not share coordinates and that mapping between the 
two types of spacetimes is vital for e.g. obtaining mock 
N-body velocity fields. 

Tildes are used to indicate fiducial spacetime points. 


Newtonian N-body simulations are based on FLRW 
backgrounds. In the case of the swiss cheese models 
studied here, the background model is the EdS model 
so this will be used as the underlying FLRW model for 
the mock N-body data. Considerations regarding the 
choice of background for the onion models are discussed 
in section VI and appendix C. 

With the background model being the EdS model, the 
mock N-body data is constructed by using the following 
relations between the EdS and Szekeres coordinates: 

t . tsz ^eds 
P Psz — Peds 
Q •— Qsz — Qeds ( 6 ) 

^^eds 

^^sz'\/ llrr,sz / ^^eds y/ 9rr,eds 

Jo 

These four equations make up a map between the 
Szekeres and EdS/N-body spacetimes. is the 

rr-component of the metric tensor in stereographic 
coordinates. The subscripts ”eds” and ”sz” are used to 
denote EdS and Szekeres coordinates respectively. The 
four equations making up the map are coupled since 
ffrr.sz depends on all four coordinates. 

Note that p and q are related to spherical coordinates 
differently in the Szekeres and FLRW spacetimes (see 
e.g. [76]). Hence, even though , fe) = (Peds,geds), it 
is generally not true that {9^^, = (^eds, <^eds)- 



Equation (6) is used to map the Szekeres density 
field onto the EdS N-body background. Aside from 
the density field, the corresponding potential, ip, will 
also be needed as it enters into the perturbed FLRW 
metric in the Newtonian gauge which is the basis for 
the ray tracing scheme This can be obtained from 


® In standard perturbation theory and thus standard ray tracing 
schemes, the coordinates of the perturbed FLRW metric are sim¬ 
ply set equal to those of the FLRW background. Thus no map is 
needed for going between these two spacetimes/coordinate sets. 


the mapped density field through the Poisson equation 
y 2 ^ _ 47 rGa with the overdensity, 5p, defined as the 
difference between the inhomogeneous ” N-body” density 
and the background density, i.e. 5p := pNbody — Pbg- 
The velocity field is needed in order to compute the 
observable redshift and the Doppler convergence (see 
subsection IV B 1 for a definition of the latter). The 
angular velocities are negligible and only the radial 
velocity, Vr, needs to be computed. The radial velocity 
can e.g. be computed from the following equation 
obtained by taking the time derivative of the r-equation 
in the map given above: 


y 

dt 


dpr{i,rsz, 


P, q) 


— (a(t)reds) = a,t(0f^eds + a{t)Vr{feds) 

( 7 ) 


dpr denotes the proper radial distance from the origin 
and the left-hand side of this equation is to be evaluated 
in the Szekeres spacetime. Vr on the right-hand side is 
the radial velocity field of the mock N-body data at the 
EdS spacetime position {t,feds,P,q) corresponding to 
the Szekeres spacetime position (t,fsz,p, g). 

Using Swiss cheese models implies involving one more 
coordinate system than the two discussed above. In 
particular, the standard ray tracing scheme is partially 
based on the global coordinate system of the swiss cheese 
model. On the other hand, the individual structures 
in a swiss cheese model are described using the local 
coordinate systems of the individual Szekeres single void 
model spacetimes making up the swiss cheese model. 
Note that the mapping described above, in principle, 
pertains to global swiss cheese coordinate systems. The 
difference between using global and local coordinate 
systems for the mapping is however negligible. 

In the following, global coordinates will be indicated 
by a subscripted ”g”. Coordinates without such a 
subscript are local coordinates. For onion models, local 
and global coordinates are the same. 


IV. LIGHT PROPAGATION AND WEAK 
LENSING 


Light rays move along null geodesics which can be de¬ 
termined by solving the geodesic equations obtained from 
the metric. In order to study weak lensing, the defor¬ 
mation of images along null-geodesics, also bundles of 
light rays must be studied. This can be done by solv¬ 
ing the screen space projected geodesic deviation equa¬ 
tion, where screen space is the two-dimensional Euclidean 
space orthogonal to the ray direction as seen by the ob¬ 
server. The Sachs formalism (based on [79]) gives a con¬ 
venient form of the formal solution through the Weyl and 
Ricci tensors. This formalism is widely used and detailed 
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descriptions can be found in e.g. [ 8 , 80] . Here, only a 

brief description of the Sachs formalism and the resulting 
ray tracing schemes are given. 


A. The Sachs formalism 


The screen space deviation vector, describes the de¬ 
viation between geodesics in a bundle of light rays. The 
Latin indices a, 6 , c G {1,2} denote screen space compo¬ 
nents while Greek indices a,/3,/x, i/,... G (0,1,2,3} will 
be used to denote spacetime components. The devia¬ 
tion vector is related to its initial derivative through a 
Jacobi relation i.e. ^(A) = D{X)^o, where dots denote 
differentiation with respect to the affine parameter, a 
subscript ” 0 ” is used to denote present time, and the 
screen space matrix D represents the Jacobi map. D is 
the central quantity for describing deformation of the ray 
bundle image; yj\ det(£))| gives the relation between the 
areas described by and Setting the initial con¬ 

ditions such that A is initially the local proper distance 
from the observer, becomes the solid angle element 
and \/| det(iA)| becomes the angular diameter distance 
Da along the fiducial ray of the ray bundle (see also e.g. 
[80]). Further more, if the cosmological model of interest 
can be considered as having an FLRW background, D 
can be related to the distortion/magnification matrix A 
by the relation D — DA.bgA [52, 82], where DA.hg is the 
background FLRW angular diameter distance. Neglect¬ 
ing rotation (see e.g. [ 8 ]), the entries of the distortion 
matrix can be decomposed into a convergence n and two 
shear components, 71 and 72 : 


A = 


- K - 71 
V 72 


72 

1 - K -I- 71 


( 8 ) 


The convergence can thus be computed from the diagonal 
component of D as follows: 


-Dll + D22 

‘2DA.bg 


(9) 


The FLRW background models used in this work contain 
non-relativistic matter and a possible non-vanishing cur¬ 
vature. In such cases, the background angular diameter 
distance can be computed using the Mattig relation [83] : 


The indices are arranged in accordance with the Einstein 
summation convention which will also be used in the fol¬ 
lowing. 

The optical tidal matrix T can be written using R := 
— YR-fivk^k'^ and F := —\Cai}g,v{e*)°‘k^{e*)^k''^ where 
is the Ricci tensor, Capg,i, the Weyl tensor and 
:= El^ + with iff, iff spanning screen space: 


/R-iie(F) im(F) A 
Im{F) R + Re{F)) 


( 12 ) 


It is convenient to note that k^ and are null- 
geodesics orthogonal to each other; this implies that the 
Ricci-part of the Weyl tensor drops out in F so that 
F = —^TZaPfj,iy(.^*)°‘k^{e*)^k‘' (this simplihcation is also 
noted in a footnote in [80]). 

In order to obtain k along null-geodesics in Szek- 
eres models, the transport equation must be solved 
simultaneously with solving the geodesic equations 
and the parallel propagation equations for iff. These 
equations are shown in appendix A together with the 
explicit expressions for R and F for the Szekeres metric. 
The resulting set of ODEs is the same as that in [52] 
though some of the components are here given in a 
slightly simpler form e.g. because the Riemann ten¬ 
sor is used instead of the Weyl tensor when computing F. 


The convergence and the reduced shear := 
are highly important cosmological observables as they 
affect redshift-distance relations and can be used to 
trace matter (see e.g. [74, 75, 84-93] for examples of 
using K and 7 in relation to observations). Especially 
with regard to the onion models it is worth mentioning 
that the quantities k and 7 are traditionally defined 
using linear perturbation theory, and their relations 
to observables is conditional to this regime. Equally 
important is it that the relation between A and D 
requires a background so that Da, hg is defined. 

In settings where the linear perturbation scheme is ill- 
suited, the algebraic relations between the components 
of D can still be used to define quantities generalizing 
K, 7 i and 72 . These quantities’ relation to observables is 
not a priori as in the linear limit (see also [52]). In this 
work, the terms ’’convergence” and ’’shear” will refer to 
the generalized definitions of k and 7 without restricting 
to linear perturbation theory. 


DA.bg = 


2 c 


RoD^.o(1 + ^)^ 

^m,0^ (2 ^-y/l J- qZ 1 


( 10 ) 


The components of D along a fiducial ray in a ray bundle 
can be obtained by using the transport equation: 


D^ = T^Dl 


( 11 ) 


^ In connection with [80], see also [81]. 


B. Standard ray tracing 

In this section, the standard ray tracing expression for 
K is introduced. For more details, see e.g. [74, 84]. 

Ray tracing combines the perturbed FLRW metric 
in the Newtonian gauge with the non-linear density 
contrasts and velocity fields obtained from N-body 
simulations. The standard recipe for going between 
Newtonian N-body results and a relativistic description 
of these is thus in line with the simple dictionary given 
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in [94] which gives a recipe for going between Newtonian 
N-body data and a relativistic spacetime with the line 
element: 


= -c^(l + 2tl))dt^+ 

('ll! 

0^(1 — 2 '(/') {dr^ + r^dO^ + s\T?{9)d(f^) 

The scale factor is normalized to 1 at present time, 
i.e. a{to) = 1. Following the recipe, the dimensionless 
gravitational potential is obtained through the Poisson 
equation, 6p, using the non-linear over¬ 

density of the given N-body data. The velocity field 
is needed e.g for computing the observable redshift. It 
is given by cx (l,u*), with u* the spatial velocity 
field of the N-body data. In the case of the Szekeres 
models studied here, the angular velocities are negligible 
and Ufi « ^(—c^(l + 2ip),a'^{l — 2tlj)vr,0,0), with 

V = ^c^{l + 2i;)-a^v^{l-2'ilj). 

The metric corresponding to the line element in 
equation (13) combined with the non-linear velocity and 
density fields will in the following be referred to as the 
Newtonian gauge metric. 


R and F can be computed and used to obtain ODEs for 
the components of 74. It is, however, standard instead 
to write the equations in terms of the deformation 
matrix A = —D and use Green’s method to obtain 

r-'A.bg 

implicit integral equations for its components. The 
resulting equations can be combined with equation (9) 
to obtain the following expression for the gravitational 
convergence: 


KS 


47rG 


dr' a^5p 


{rg - r'„)r'„ 


(14) 


This integral expression is obtained by neglecting lens- 
lens coupling (see e.g. [95]) and by using the Limber ap¬ 
proximation (see e.g. [68, 96]). The integrals are along 
the path of a given light ray and are usually solved by in¬ 
troducing the Born approximation, where rays are traced 
as radial null-geodesics in the background model. 

The convergence gives the most important contribu¬ 
tion to alterations of redshift-distance relations in N- 
body simulation models where it is standard to use the 
approximation « Da, hgi^ — ks). (1 —k^) is the lowest 

order approximation to the determinant of D since the 
shear is usually several orders of magnitude smaller than 

KS. 


1. Doppler convergence 

The total first order convergence consists of several 
other contributions than just the gravitational con¬ 
vergence. These other contributions e.g. include the 
well-known Sachs-Wolfe and integrated Sachs-Wolfe 
effects which are generally subdominant to the gravita¬ 
tional convergence. There is also a Doppler contribution. 


Ky, which arises because peculiar velocities affect the 
redshift along light rays. This Doppler contribution to 
the convergence and redshift-distance relations has been 
known for several decades (see e.g. [97] for an early 
derivation), and has been included in several more recent 
studies {e.g. [72, 98-104]). The Doppler convergence is 
typically not included when ray tracing through N-body 
simulations. This is presumably because the effect is 
important mostly at small redshifts, and because the 
effect should become negligible when averaged over an 
infinite number of observations in a universe satisfying 
the Copernican principle. The Doppler contribution will 
be included here since it is automatically included in the 
geodesic deviation equation and thus will be portrayed 
in the exact convergence. 

The Doppler convergence can be computed as (see 
[98]): 

KyC = [1- ^ -- ) n-(v-Vo)-kn-Vo (15) 

V aAVo-V)J 

The vector n is the photon direction vector in the source 
frame i.e. it is the spatial part of the unit 4-vector 

oc evaluated along the null geodesics, rj 

denotes conformal time. 

The velocity field v is the spatial part of the normal¬ 
ized 4-velocity . The sign of v is determined from the 
global coordinate system and can in swiss cheese models 
differ from the local single void sign of . 

Note that in equation (15), the scale factor etc. 
are evaluated at the spacetime position corresponding 
to the observed redshift, while equation (14) for the 
gravitational convergence is evaluated in accordance 
with the background redshift. 


2. Summary of the ray tracing method 

In the following two sections, a simple and popular 
approximate ray tracing method is compared with exact 
light propagation and with an extended ray tracing 
method which includes the Doppler convergence. For 
convenience, the approximations included in the ray 
tracing methods are summarized here. 

The convergence is computed by adding the lowest 
order gravitational convergence, ks, to the Doppler 
convergence, The gravitational convergence is 

computed using the Born and Limber approximations 
and by neglecting lens-lens coupling. In standard ray 
tracing schemes, the Doppler convergence is neglected. 
Thus, in standard ray tracing schemes, the angular 
diameter distance is computed as Da = 74^ bg(l — ks) 
while in the extended ray tracing scheme included here, 
it is computed as Da = 74^_bg(l — ks — Ky). 

The convergence will be portrayed as a function of 








the observable redshift given by Zobs + 1 = ■ The 

subscript ”e” denotes time of emission. 

The potential ip turns out to be insignificant for 
the results. This is only the case since ks is computed 
using the Limber approximation and thus is computed 
independently of ip. The only place ip plays a role is 
thus when lowering indices in the computations of Ky 
and Zobs- The most significant ’’first order” quantity 
modifying these from their background values is the 
velocity field, not ip. By using a few of the studied LTB 
models, it has been tested that simply neglecting ip 
altogether is of no significance for the results obtained 
here. The results presented below have thus all been 
obtained under the approximation of setting ip = 0. 

V. SWISS CHEESE MODEL RESULTS 

In order to compare results obtained using the exact 
formalism and the approximate ray tracing schemes, 
these methods will be used to study the convergence 
along individual rays in swiss cheese and onion models. 
A direct comparison between the exact and approximate 
results requires that the exact and ray traced rays 
are initialized equivalently. Since the exact ODEs are 
solved backwards in time, this implies that the observer 
position and line of sight in the two spacetimes must be 
equivalent according to the map in equation (6). 

When the purpose is studying light propagation, 
Swiss cheese models are most simply constructed on 
the fly. In such construction schemes, a ray will be 
moved from one single void model to another once it has 
reached a specific coordinate distance from the given 
local structure (see also [8]). In general, the exact and 
ray traced rays do not follow the same paths. Since the 
Swiss cheese models are made on the fly, this implies 
that the swiss cheese models that the exact and ray 
traced rays propagate through are unrelated. Clearly, 
this leads to a rather contrived comparison between the 
exact and ray tracing results. To avoid this, when a ray 
traced ray is moved to a new structure, its position and 
direction will always be chosen to be the initial position 
and direction that the corresponding exact ray had 
when it reached that specific structure ® . This way, the 
exact and ray traced rays will move through the ’’same” 
swiss cheese model, but the rays will generally traverse 
each structure along different paths - depending on how 
good an approximation the Born approximation is. The 
swiss cheese results are thus to be considered indicative 


® Within a given redshift interval, the ray traced ray may traverse 
more structures than the exact ray. When traversing such final 
structures, the ray traced ray is moved to a ’’new” structure 
simply by turning it on the spot and tracing it radially back 
towards the structure it came from. 


of accumulative effects of local inhomogeneities. 

The standard ray tracing scheme described in section 
IV B is adapted to radial rays. In practice, ray tracing 
through Swiss cheese models thus becomes much simpler 
if rays are always initialized as radial at the beginning 
of a new single void structure. The ray tracing results 
shown below have therefore been obtained by following 
the procedure of moving a ray to a new structure by 
turning it around and sending it radially back towards 
the structure it just left. Especially the part of ini¬ 
tializing rays as radial at each ’’new” structure clearly 
makes the paths very special. The consequent bias of 
the results induced by studying only special rays implies 
that it is inappropriate to compare with background 
results. This does not affect the integrity of the study 
as the point is not to compare with background results 
anyway. The point here is to study the precision of the 
standard ray tracing scheme with the added Doppler 
convergence as well as to study implications of neglecting 
the Doppler convergence. Thus, for the results obtained 
from the swiss cheese models based on non-symmetric 
Szekeres models, the special paths do not compromise 
the study much; if the ray traced rays cannot reproduce 
exact rays that are initialized as radial then surely they 
cannot reproduce more general rays that are initialized 
with arbitrary impact parameters at each structure. 

The consequence of only studying radial rays is much 
more significant in the case of the spherically symmetric 
LTB models; in these models, rays initialized as radial 
will stay so. This implies that the exact and ray traced 
rays will automatically move along the same spatial 
path. In order to see how (dis-)similar the exact and 
ray traced paths along non-radial rays in LTB models 
are, the Sachs formalism has been used to obtain an 
approximation scheme for computing the convergence 
along such non-radial rays. The scheme employs the 
Born approximation and should yield results that are 
no more precise than the regular ray tracing scheme 
described in section IV B. The details of the scheme are 
explained in appendix B. 

Unless otherwise stated, the observer is placed at 
the center of a void. The results are always plotted 
against the observable redshift and light rays are prop¬ 
agated/traced through swiss cheese models until they 
reach a redshift of Zobs = 0.6. 


A. LTB results 

In this section, the results from studying light prop¬ 
agation in a swiss cheese model based on LTB single 
void models are shown. As discussed above, the results 
shown here pertain to rays that propagate radially 
through consecutive LTB single void structures. This 
makes the rays very special since is guarantees that the 
exact and ray traced rays move along the same spatial 
paths. Non-radial rays in LTB models will be discussed 
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in subsection V A 1. 


As seen in figure 3(a), the local density distributions 
along the approximate and exact geodesics are the same. 
This indicates that the combined (mock) N-body data 
and Newtonian gauge metric yield a good approximation 
of the exact observable redshift. Note that it is impor¬ 
tant to use the approximate observable redshift; simply 
using the background redshift leads to clear deviations 
between the exact and ray tracing results. 

The convergence along the ray is shown in figure 3(b). 
There is an overall good agreement between the exact 
convergence and the sum of the Doppler and gravita¬ 
tional convergences. The figure also shows that the 
Doppler convergence clearly dominates the convergence 
at low redshifts while the gravitational convergence 
dominates at higher redshifts. The former is in agree¬ 
ment with the results obtained in [72] where the Doppler 
and gravitational contributions to the convergence were 
studied in Szekeres single void models. In [72], mock 
N-body data was not employed when computing the two 
contributions to the convergence. Instead, regular linear 
perturbation theory was used. It was found that, at late 
times when the Szekeres structures had gone non-linear, 
linear perturbation theory lead to an underestimation 
of the Doppler convergence by approximately 20%. 
Here, it is found that this underestimation vanishes 
when linear perturbation theory is combined with 
non-linear Newtonian N-body velocity fields. This is 
emphasized in figure 4 which includes close-ups along 
the convergence curves. Not until approximately z > 0.5 
do the exact and ray tracing convergences begin to 
differ notably from each other - and even then, only by 
a few percent. The difference between the exact and 
approximate convergences grows slightly with ri, and m 
(defined in section II A). It is virtually nonexistent in 
the Tb = 30Mpc cases, and for the m = 2 cases also 
when rb is larger. The figures shown here are for LTB 
models with m = 4 and rb = 60Mpc. 

The difference between the exact and ray traced 
convergences are at least partially numerical artifacts. 

The distinct oscillations of the Doppler convergence 
appear as a manifestation of the Doppler convergence 
being a measure of the local velocity field. The grav¬ 
itational convergence is an accumulated effect of the 
matter distribution along the given geodesic. The 
gravitational convergence is thus much less affected by 
individual local variations than the Doppler convergence 
is, though close-ups (not shown) of the gravitational 
convergence do reveal small ripples accounting for the 
local fluctuations in the density field. 


It is clear from figure 3(b) that the Doppler con¬ 
vergence, in principle, will affect observations. As 
a simple study of the significance of these effects, 
the different convergences have been used to com¬ 


pute the luminosity distance as a function of redshift 
which has then been used to obtain estimates of the 
reduced Hubble parameter, /i, and the deceleration 
parameter, go- The luminosity distance is computed 
as Dl = Da{ 1 + zY = Da, hg,Y ~ '^)(1 + In 
low redshift interval of 0.01 < z < 0.06 ® , 200 points 
along the redshift-distance relations are fitted to the 
low-redshift FLRW approximation Dl{z) « j|-z with 
h as the fitting parameter. In the remaining redshift 
interval of 0.06 < z < 0.6, 1000 points are fitted to the 
FLRW approximation Dl{z) « ^z (l-I-I-^z) with 
go as the fitting parameter, and using the value of Hq 
obtained from the low-redshift fit. 

To illustrate the significance of void size, the fit¬ 
ting results are shown for two void sizes, namely 
Tb = 30Mpc, lOOMpc. The fitting was done with pyx- 
plot^ and the best-fit values together with 1 standard 
deviation intervals are shown together with the plots in 
figure 5. 

There is a distinct difference in the estimates of h 
based on the ks fits compared to the fits based on Kexact ■ 
This again affects the estimates of go. The difference in 
the two estimates of h is only approximately 0.1%, which 
is quite small compared to uncertainties of observational 
results. For instance, in [106] the Hubble constant is 
estimated to be Hq = 73.8 ± 2.4km/s/Mpc i.e. the 
uncertainty is 3.3%, and the authors express no hope 
of their procedure in the future to yield results with a 
precision below approximately 1%. 

The fits obtained using the sum ks + Ky agree well 
with the exact results. In particular, the two best-fit 
values of h agree within approximately 0.02% for the 
rb = 30Mpc model and 0.03% for the rj, = lOOMpc 
model. The void size clearly affects the determination of 
both and go.exact, go,Ki-i-K„ both regarding 

best-fit values and their uncertainties. 

A notable feature of the fits is that the uncertainties 
are markedly higher for the fits obtained using Kexact and 
the sum Ky + ks than for those obtained using just ks 
and background values of Dj^. This finding agrees with 
what should be expected from figure 3(b); the oscillating 
Doppler convergence leads to larger uncertainties in the 
low redshift estimates and to deviations of the best-fit 
values from those obtained from the background Dj^ 
and Ks- At higher redshifts, the Doppler convergence 
becomes less important and the uncertainties in the 
best-fit values become more similar. 

As is seen in especially figure 5(b), the imprints of the 
Doppler convergence oscillations on observables can be 
very prominent. In fact, the Doppler convergence is so 
prominent at low redshifts that it in [99] was studied to 
what extent measures of the Doppler convergence can 
be obtained observationally with the prospect of using 


® This interval is chosen based on figure 1 of [105]. 
^ http://pyxplot.org.uk/ 
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Density convergence 




(a) (b) 

FIG. 3. Density and convergence along a light ray in a swiss cheese model based on compensated LTB single void models with 
m — A and Vb = 60Mpc. The ray path has been computed using both the Born approximation and the geodesic equations of 
the exact LTB metric. The two resulting density distributions are indistinguishable in the figure. Four convergences are shown. 
The line labeled shows the sum of the Doppler convergence and the gravitational convergence obtained using the 

simple ray tracing scheme. The two components are also shown individually. The line labeled ’’Wexact” is the exact convergence 
as obtained from the exact set of ODEs. The figure includes a close-up of the region where the gravitational convergence begins 
to be non-negligible. 


convergence 



0 0.2 0.4 0.6 

^obs 



O.fi 


FIG. 4. Gonvergence along a radial ray in the swiss cheese model based on the LTB single void model specified by m = 4 and 
Vb = 60Mpc. The line labeled -(- ns" shows the sum of the Doppler convergence and the gravitational convergence obtained 
using the simple ray tracing scheme. The line labeled ’’Kexact” is the exact convergence as obtained from the exact set of ODEs. 
The figure includes close-ups at different redshift intervals. 


it to e.g. study peculiar velocity fields. 


1. Non-radial rays 

It was shown above that the simple ray tracing scheme 
with the convergence computed as -I- ks yields an 
almost exact reproduction of the exact results. This 
exact reproduction does not depend on the rays being 
radial. To show this, the Born approximated Sachs 
formalism described in appendix B has been used to 
compute an approximate convergence along non-radial 


rays in different LTB models. Results are shown in figure 
6 for the LTB single void model specified by m = 4 and 
Tb = 60Mpc. As seen, the exact paths are exactly traced 
out by the Born paths. Aside from a minor oscillation 
at Zobs ~ 0, the convergences are also reproduced very 
precisely. These oscillations are due to changes in the 
velocity field which have, for simplicity, been neglected 
in the ray tracing method described in appendix B. The 
oscillations are, however, described very well by the 
standard ray tracing method described in section IV B. 

Before moving on to study non-symmetric Szekeres 
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Luminosity distance, = SOAIpc 


Luminosity distance, = lOOMpc 




(b) 


Luminosity distance, — SOAIpc 


Luminosity distance, rf, = lOOAIpc 




(c) 


(d) 


FIG. 5. Luminosity distance plotted as a function of the observed redshift and fitted to low-redshift FLRW approximations 
of Dl{z). The results from four luminosity distances are shown corresponding to the exact convergence, the sum of the 
gravitational and Doppler convergences, the gravitational convergence and the background solution. The results are shown for 
LTB Swiss cheese models specified by m = 4 and both rt = 30Mpc and rb = lOOMpc. 



- Born, /cq — — ^ 

" exact, ^0 — ^ ^ 

. Born, ^0 — — f 

^ exact, ^0 — — ^ 

- Born, ^0 ~ “c 

■ exact, ^0 — ~ I 


FIG. 6. Density and convergence along non-radial light rays in an LTB single void model specified by m = 4 and rt = 
60Mpc. The observer has been placed at {t,r,9,(j>) = (to, 70Mpc, S/Stt, S/Stt), looking in various directions determined by 
{k^ ,kF ,k'^) = (—1, fco, fcg, 0), where k^ is varied as fcg = and feo is determined from the null-condition. Note 

that fcg does not have the same units as 1 so it is only the numerical value of c that is implied in the values of fep. Since the 
LTB model reduces exactly to the EdS model outside the structure, these initial conditions are for both the exact Szekeres 
coordinate system and the ray tracing coordinate system. The exact and ray traced rays are labeled by ’’exact” and ’’Born” 
respectively. 


models, it should be pointed out that the results of this section show that the LTB models can generally 
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FIG. 7. Density along geodesics in swiss cheese model based on compensated non-symmetric Szekeres single void models. The 
geodesics have been computed using both the Born approximation and the exact geodesic equations according to the exact 
Szekeres metric. 



be described extremely well using the standard ’’pertur¬ 
bation” schemes used for studying N-body data. The 
quotation marks around ’’perturbation” emphasize that 
the ray tracing scheme is not strictly based on linear 
perturbation theory as the velocity and density fields 
are non-linear. This result is in agreement with the 
results of [76] where it was shown that the Newtonian 
gauge metric could reproduce light propagation in a 
specific uncompensated LTB single void model, while 
it was insufficient for describing light propagation in 
non-symmetric models. 


B. Non-symmetric Szekeres results 

In this section, results from studying exact and 
ray traced rays in swiss cheese models based on the 
non-symmetric Szekeres models with rf, = 60Mpc and 
TO = 2 , 6 are shown. 

Because of the anisotropy of the voids, the exact 
rays are not radial. As seen in figure 7 this results 
in clearly different densities along the exact and ray 
traced rays, with the difference being most significant 
for the TO = 2 case. In particular, the exact rays seem 
to completely avoid the central most devoid parts of the 
voids - even though they are initialized as radial at each 
single void structure. On the other hand, the ray traced 
rays are forced through the centers of each void. 

In accordance with this clear difference of local paths, 
the convergence along the exact and ray traced rays 
are also different. The convergences along the rays are 
shown in figure 8 . Aside from being different locally, the 
exact and ray traced convergences do not even converge 
towards each other at high redshifts. This is because the 
gravitational convergence dominates at high redshifts. 
Especially in the to = 2 case, the very different densities 


along the exact and ray traced paths apparently lead 
to gravitational convergences that are different. In the 
TO = 6 case, the exact and ray traced convergences have 
almost converged to each other at Zobs ~ 0.6. 

It is also notable that, unlike in the LTB case, 
the Doppler convergence does not oscillate about the 
gravitational convergence. This is e.g. seen by noting 
that the gravitational convergence is approximately 
zero at very low redshifts. If the exact convergence 
oscillated about the gravitational convergence, it would 
thus oscillate about zero at these lowest redshifts. This 
is clearly not the case. 


As an illustration of the importance of the differences 
in the exact and ray traced paths, the luminosity 
distances corresponding to the exact and ray tracing 
convergences are used to estimate h and go- The esti¬ 
mates are obtained using the fitting scheme introduced 
in section V A. The results are shown in figure 9. 

The parameter values determined from the exact and 
ray traced rays are clearly in less agreement than in the 
LTB case. Specifically, the determinations of h based 
on Kexact and Ky + Ks differ by approximately 2.9% and 
3.9% in the to = 2 and to = 4 cases respectively. The 
deviation between the exact and ray tracing estimates 
of h depends significantly on the line of sight of the 
observer. The results shown here are for a direction 
where the difference between the approximations is par¬ 
ticularly large. In many other directions the deviation is 
sub-percent. It should in this connection be emphasized 
that the point with the estimates made here is not to 
determine whether or not local inhomogeneities may af¬ 
fect determinations of h (or go) in the real universe; this 
would require averaging over many rays with variable 
impact parameters at each void. This has, for instance, 
been done in [ 8 ] and as mentioned in the introduction. 
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convergence, = 60A/pc, m = 2 


convergence, — QOMpc, m — Q 




FIG. 8. Convergence along rays in swiss cheese models based on non-symmetric Szekeres single void models. Four convergences 
are shown in each figure. The line labeled ” Kv-\-ks^^ shows the sum of the Doppler convergence and the gravitational convergence. 
The two components are also shown individually. The line labeled "/^exact” is the exact convergence as obtained from the exact 
set of ODEs. 


Luminosity distance, = QOMpc, m — 2 


Luminosity distance, — QOMpc, m — Q 




Luminosity distance, — QOAIpc, m — 2 


Luminosity distance, = 60Afj)c, m = 6 




(c) 


(d) 


FIG. 9. Luminosity distance plotted as a function of the observed redshift and fitted to low-redshift FLRW approximations of 
Dl{z). The results from four luminosity distances are shown corresponding to the luminosity distance obtained from the exact 
equations, the sum of the gravitational and Doppler convergences, the exact convergence and the background solution. 
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the conclusion of that and other studies is that when 
averaging over many rays, the distance-redshift relation 
reduces to that of the average model. Thus, even though 
there may be differences in the exact and ray tracing 
estimates of h along individual rays, the average results 
should be similar. The point with this section is thus 
that the exact paths of rays are significantly affected by 
inhomogeneities and to such an extent that it clearly 
affects observables along the individual rays. It is shown 
here that this is not described correctly by simple ray 
tracing schemes. This may be important in relation to 
interpretation of observations such as shear maps and 
CMB temperature fluctuation maps, where averaging 
over large portions of the sky is not viable. 

The differences in the estimates of are also quite 
large in the example shown here, but this is mainly due 
to the differences in the estimates of h. 

The large deviation between exact and Born ap¬ 
proximated paths is highly dependent on the level of 
anisotropy. If the anisotropy is lowered to approximately 
25% of that used here, i.e. if Q defined in section 
II is multiplied by 0.25, then the exact and Born 
approximated rays are almost indistinguishable. 


VI. ONION MODEL RESULTS 

The Swiss cheese studies indicate that, in as far as the 
Born approximation is a good approximation, the simple 
ray tracing scheme studied here is very precise. The swiss 
cheese studies also indicate that the Born approximation 
generally cannot be assumed to be a good approximation. 
In the Swiss cheese studies, the global paths of the exact 
and ray traced rays were the same since a ray traced ray 
was initialized at the same position and in the same direc¬ 
tion as the corresponding exact ray at each new structure. 
The onion models are multiple structured models ab ini¬ 
tio and not only after an on-the-fly construction scheme. 
Whereas the swiss cheese models were used to study the 
local precision of the Born approximation, onion models 
can thus also be used to study the global precision of the 
Born approximation. In addition, the convergence along 
radial rays in LTB onion models will be studied in order 
to examine the validity of using the ray tracing scheme 
in models without intrinsic backgrounds. 


A. LTB results 

The convergence along light rays in LTB swiss cheese 
models was very precisely reproduced by the simple ray 
tracing scheme that included the Doppler convergence. 
It is here tested whether or not this conclusion also holds 
true for radial rays in LTB onion models. The results 
are shown in figure 10. 

The density fields along the rays indicate that the 
redshift is approximated well by the ray tracing scheme. 


An exception is in model Olitb which has the largest 
truncation radius of the studied onion models. In figure 
10 (a) it is seen that the densities along the exact and 
ray traced rays in model Oljtb diverge slightly from each 
other at the highest part of the studied redshift interval. 
Since the exact and ray traced rays automatically move 
along the same spatial path, this indicates a deviation 
between the exact and approximate redshifts. The fig¬ 
ures also reveal that the convergence is not reproduced 
quite as precisely in the LTB onion models as in the 
swiss cheese case. The deviation between the exact 
and ray tracing convergences grows with the truncation 
radius i.e. the deviation is most prominent in model 
Olitb while it almost vanishes in model 03itb- 


It may be that a better reproduction of the exact re¬ 
sults can be obtained if the ray tracing is done on another 
global background than the EdS model. The theory pre¬ 
sented in sections III and IV B is based on the assumption 
of a flat background. Modifications needed in order to 
adapt the equations to curved backgrounds is described 
in appendix C together with discussions and examples 
regarding the use of curved FLRW backgrounds. Using 
these modified equations, the ray tracing through model 
Olitb has been done on different backgrounds. The re¬ 
sults are shown in figure 11. As seen, the accuracy of the 
reproduction of the convergence depends on the back¬ 
ground. The best reproduction is obtained when using a 
background defined by a present time density parameter 
0 := ~ 0.83. The average density along the ray 

corresponds to ~ 0.8. Using a background with a 
small negative curvature clearly yields a better reproduc¬ 
tion of the exact results although the reproduction is still 
not quite as good as along rays in swiss cheese and single 
void models. It seems likely that this is simply a result 
of the density fluctuations in the onion model oscillating 
about a ’’mean” that changes in the radial direction. In 
particular, the mean starts below, but converges towards, 
the critical density. This feature is not particularly real¬ 
istic and it seems likely that the reproduction of the con¬ 
vergence and redshift would be better in more realistic 
models or if the studies of the models used here were re¬ 
stricted to redshift intervals where the oscillations occur 
approximately about the same mean. For instance, fig¬ 
ure II shows that at redshifts below approximately 0.15, 
the change in this mean is very small. In accordance 
with this, the convergence in the same redshift interval is 
well described by ray tracing with a constant background 
corresponding to = 0.8. 

B. Non-symmetric Szekeres results 

The convergence is not correctly reproduced along 
the radial rays in the LTB onion models when an EdS 
background is used for the ray tracing. Instead, a curved 
background must be introduced. The background 
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Density along radial ray, model 01^6 


Density along radial ray, model 02itb 


Density along radial ray, model O'Sub 




(d) 


(e) 


(f) 


FIG. 10. Density and convergence along exact and ray traced radial rays in the onion models Olitb, 02itb and OSitb. The exact 
convergence is shown together with the sum of the Doppler and gravitational convergences. The gravitational convergence 
is also shown alone but its deviation from zero is too small to be visible in the figure. For model 03itb, a close-up of the 
convergence has been included to highlight that the exact and approximate convergences are not exactly identical. 


necessary for obtaining accurate reproductions depends 
on the density fluctuations along the given ray. Even for 
a central observer in the nonsymmetric Szekeres models, 
these density fluctuations depend on the line of sight. It 
is thus not entirely clear which background to choose as 
the global background for the ray tracing. In addition, 
as discussed in section VIA, the introduction of the 
damping of the density fluctuations makes the models 
unsuitable for studies where FLRW backgrounds are 
introduced. However, a somewhat precise reproduction 
of the exact results is obtained along the entire rays 
when ray tracing on an EdS background in models 02itb 
and 03itb- Thus, as an initial consideration, the exact 
and approximate convergences along individual rays in 
models 02sz and 03sz will be compared here while model 
OIsz will not be studied. A more thorough, future study 
of the effects of nonsymmetric, interacting structures 
when ray tracing with fixed FLRW backgrounds, could 
e.g. be based on the models introduced in [108, 109]. 

The density distributions and convergences along 
individual rays initialized at an observer at the origin 
in model 03sz are shown in figure 12. The results 
are shown along the three lines of sight defined by 
(0sz.o,<^sz.o) = 3/5^(I,I), (0sz.o,</>sz.o) = (3/57r, II/IO^) 
and (0sz,o, <^sz,o) = (3/57r, 43/307r). The first direction 
is close to the direction of highest anisotropy while the 
second is approximately orthogonal to it. The last line 


of sight is in a direction in between the first two. Note 
in figure 12 that the exact and ray traced rays deviate 
less from each other the closer the line of sight is to the 
line of highest anisotropy. This is not surprising since 
the density is close to being symmetrically distributed 
about this line. 

The three sets of results have been used to esti¬ 
mate the reduced Hubble parameter, h, as in earlier 
sections. The resulting estimates of h and the corre¬ 
sponding uncertainties are shown in figure 13. The 
estimate of h is approximately 1% higher when based on 
Kexact than when using Ky + Kg. 


A single ray has been studied in model 02sz. 
The ray is initialized along the line of sight 
(^sz, 0 ) ?^>sz,o) = (3/57r, 43/307r), i.e. in between the 
lines of highest and lowest anisotropy. The convergence 
and density field along the exact and ray traced rays 
are compared in figure 14. As with the 03sz results, 
the exact and ray traced rays clearly do not move along 
equivalent spacetime paths. 

In model 02sz, the inhomogeneities in the density 
distribution are truncated at a large enough distance 
from the origin for it to be reasonable to make an esti¬ 
mate of the present time deceleration parameter, qq. h 
is estimated from 200 data points in the redshift interval 
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Density along radial ray in model Olub with = 0-65 


Density along radial ray in model Olub with flm.o = 0.8 




Density along radial ray in model Olub with o = 0.83 


Density along radial ray in model Olub with = 0.95 



FIG. 11. Density and convergence along exact and ray traced radial rays in the model Olitb- The ray tracing has been done on 
different backgrounds as indicated in the individual figure titles. The average density of model Olitb along the ray corresponds 
approximately to r2m,o = 0.8. 
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Density along ray in model OSs^ with [dg^.o, <l>9z,o) = 3/57r(l, 1) 



(a) 


Density along ray in model 03.,- with (0,2,o><^'s2,o) = {3/57r, 43/307r) 



Zobs 

(b) 


Density along ray in model 03^2 with (0s2,o>‘/’s2,o) = (3/57r, II/IOtt) 



(c) 


convergence along ray in model 03 .,2 with ( 0 , 2 , 0 .0*2, 0 ) = 3/57r(l, 1) convergence along ray in model 03s2 with (0s2,o, <i>sz,o) = (3/57r,43/307r) convergence along ray in model 03*2 with (0s2.o, ^s 2 .o) = (3/57r, II/IOtt) 





FIG. 12. Density and convergence along exact and ray traced rays in onion model 03sz. Four convergences are shown: Doppler, 
gravitational, sum of these two, and the exact. The gravitational convergence is negligible in the low-redshift interval studied 
here, so the lines showing the Doppler convergence and the sum of the Doppler convergence and the gravitational convergence 
are indistinguishable. 


Luminosity distance, model 03^2 with (0s2,oi‘/’s 2 ,o) = 3/57r(l,l) Luminosity distance, model 03,2 with (0s2,o> ^ 22 , 0 ) = (3/57r,43/307r) 


Luminosity distance, model 03.,- with (0.,2,o><^',2,o) = {3/57r, II/IOtt) 





(a) 


(b) 


(c) 


FIG. 13. Luminosity distance plotted as a function of the observed redshift and fitted to the low-redshift FLRW approximation 
for Dl{z). The results from four luminosity distances are shown corresponding to the exact convergence, the sum of the 
gravitational and Doppler convergences, the gravitational convergence and the background solution. The results are shown 
along three different lines of sight for a central observer in onion model 03az- 


0.01 < 0 < 0.06 while qo is estimated from 200 points 
along the redshift interval 0.06 < z < 0.12. Here, the 
interval and number of data points used to estimate Qq is 
much less than in the swiss cheese cases which explains 
the larger uncertainties in the results shown in figure 15. 
The estimates of h based on /texact and the sum Hy + ks 
are quite similar with a difference of approximately 
0.3%. The difference between the two corresponding 
estimates of go are significantly different; the estimate 
based on Kexact is more than twice the estimate based on 
the sum Ky + ks- This result is very different from that 
obtained by studying the swiss cheese models where the 
different estimates of go were primarily due to significant 
differences in the estimates of h. The difference in the 


estimates obtained here is in line with the convergences 
along the two rays shown in figure 14 where Ky -|- ks is 
seen to be more negative than Kexact- The difference 
in the go estimates may merely be a manifestation of 
the differences in the density distributions along the 
exact and ray traced rays. Considering the results found 
in section VIA it seems likely that the differences are 
partially due to the apparent insufficiency of the ray 
tracing scheme to describe light propagation in onion 
models. 


Lastly it is noted that light propagation through 
onion models has its limits regarding how realistic it is; 
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Density along ray in model 02„ with (0„,o, </>.*,o) = (S/Srr, 43/307r) convergence along ray in model 02„ with (6l„,o, 0„,o) = (3/57r, 43/307r) 




(b) 


FIG. 14. Density and convergence along exact and ray traced rays in onion model 02sz. Note that tci « 0 so that the lines 
representing and the sum + ks are indistinguishable in the figure. 


Luminosity distance, model 02^2 with {Osz, 0 jfpsz,o) = (3/57r,43/307r) Luminosity distance, model 02sz with 0s2,o) = (3/57r, 43/307r) 






(a) 


(b) 


FIG. 15. Luminosity distance plotted as a function of the observed redshift and fitted to low-redshift FLRW approximations 
for Dl{z). The results from four luminosity distances are shown corresponding to the exact convergence, the sum of the 
gravitational and Doppler convergences, the gravitational convergence and the background solution. 


as was seen in figure 7, real light rays will avoid the most 
underdense regions of a void. In onion models, the light 
rays are less free to bend around the most underdense 
regions since these are spread out over large angular 
regions. 


VII. CONCLUSIONS 

By studying both radial and non-radial rays in LTB 
Swiss cheese and single void models it was shown that 
the sum of the gravitational and Doppler convergences 
reproduces the exact convergence very precisely. The 
Doppler convergence is not included in standard ray 
tracing schemes, but in accordance with other stud¬ 
ies, it is here found that the Doppler convergence is 
a very significant contribution at low redshifts. Its 
significance for parameter determinations was here 
illustrated through estimates of the reduced Hubble and 
deceleration parameters. It was found that the Doppler 
convergence has a distinct yet rather small effect on 


especially the uncertainty of these estimates. 

In the non-symmetric Szekeres models, exact rays 
and ray traced rays do not propagate through equivalent 
portions of spacetime. The difference between exact and 
ray traced paths is very dependent on the void profile. 
The more square-like the void profiles are, the closer the 
Born approximation comes to being exact. In the case 
where the void profiles are very cone-like, the densities 
along the ray traced rays and the exact rays are very 
different because the exact rays avoid the central part 
of the void. Even in this case, the consequence for e.g. 
determinations of the reduced Hubble and deceleration 
parameters generally seems to be limited. However, 
along some particular lines of sight, the differences in the 
estimates of h based on the exact Da and that obtained 
using Ky + Ks can become ^ 3 — 4%. The effect might 
be more significant for other quantities. It would, for 
instance, be interesting to see if there are any significant 
differences between CMB temperature fluctuation maps 
obtained using exact rays and Born approximated rays. 
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Many weak leasing studies employ more sophisticated 
ray tracing schemes than the one used here. In par¬ 
ticular, the Born approximation is often replaced by 
multiple-lensing techniques where the ray traced rays 
are bent at a discrete number of leasing planes. It would 
be interesting to see to what extent multiple-lensing 
techniques increase the agreement between the exact 
and ray tracing light paths. The study in [76] indicates 
that even if the full Newtonian gauge metric is used to 
propagate rays, the resulting ray paths will generally 
be different from the corresponding exact ray paths. It 
is thus not clear that using multiple-lensing techniques 
would result in significantly more correct ray paths than 
the Born approximation. 

The structure formation in swiss cheese models is 
very special since each single void structure is inhibited 
from interacting with the other structures in the swiss 
cheese. To test if this has any significance for the 
precision of the ray tracing scheme, onion models have 
also been studied. 

In the case of LTB onion models, the accuracy and 
precision of the ray tracing results are not quite as good 
as in the case of the LTB single void and swiss cheese 
models, and for one of the studied models, obtaining 
a somewhat good reproduction requires introducing a 
curved FLRW background. The study based on radial 
rays in LTB onion models indicates that the damping 
of the density fluctuations in the models compromises 
the study. Generally though, the results are i overall 
agreement with those found for the swiss cheese mod¬ 
els. A more thorough study of effects of interacting, 
nonsymmetric, dynamical structures on ray tracing on 
fixed FLRW backgrounds may be in order and could 
for instance be done using the models introduced in 
[108, 109]. 
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is to obtain k along light bundles. This can be achieved 
by simultaneously solving the geodesic equations, the 
parallel propagation equations for the components of 
and the equations for the components of the Jacobi 
matrix D. ODEs for this were obtained in [52] and later 
used in [8]. The ODEs given below are the same as 
these, with the notation used here being slightly different. 


To simplify the notation, the metric components 
will be written using functions defined according to the 
line element written as: 


ds^ = —c^dv 


1 — fc(r) 


■dr^ 




;{dp^ + dq^) 


E(r,p,qY 

= -c^dt^ + R{t,r,p, q)dr"^ + F{t,r,p,q) {dp^ + dq^) 

(Al) 


The geodesic equations corresponding to this line 
element are: 

+ m) (A2) 

F = ^ (-2i?fc’- + RAkn^ + Fr [(F)2 + (A:«)"]) 

(A3) 

k^=l^ {-^Fk^ + RAkl^ + F,p [{kPy + (F) 2 ]) 

(A4) 

A = ^ (-2W + R )2 + F, [{kPf + (F)2]) 

(A5) 

The parallel propagated screen space basis vectors, E^, 
are needed since they enter into the ODEs for D. ODEs 
describing this parallel propagation are simply obtained 
by substituting ^ {k^F^ + k''E^) for k^^k’' in the geodesic 
equations with the resulting equations being: 

^a =iR,tF:k^ + F,t [EPkP + F«A:«]) (A6) 


Appendix A: The Sachs formalism specified to 
Szekeres models 

In this appendix, the Sachs formalism described in sec¬ 
tion IV A is specialized to the Szekeres model. The goal 


2i?F: = -i?,* [F*F + Elk*] - R,rk^E:- 
R,p [ElkP + El¥] - A,, [El¥ + El¥] + F, [F^F + F^F] 

(A7) 


htt p: / / WWW. gnu. org / software / gsl / 


2EEI = -Ft [El¥ + Elk*] + F^- 

F, [¥El + ¥El] - F^^El¥ + E^^El¥ - F, [El¥ + El¥] 

(AS) 
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Explicit expressions for R and F need to be obtained. 
The expression for R is easy to obtain; the Einstein equa- 
—tion can be written as 7?.^^ = + SttGw Com- 

[¥El + - Fp [FF« + k'^E^] - E^qk'^E^ + F^PF^bining this with the null condition k>^k^ = 0 and noting 

(^A9) that the Szekeres spacetime is given in a comoving space- 
time foliation such that = (1, 0, 0,0), this leads to: 


The last ODEs needed are those for the components of 
D which by equation (11) are seen to be: 


Ziii = (R-Fe(F))Fii + Jm(F)Fi 2 (AlO) 


R = - 47 rGArp(fc *)2 (A 14 ) 


Fi 2 = (R-Fe(F))Fi 2 + /m(F)F 22 (All) 


F 21 =/m(F)Fii + (R-tFe(F))F 2 i (A12) 

1)22 =/m(F)Fi2 + (R + Fe(F))F22 (A13) 


To obtain an explicit expression for F, the Riemann 
components are needed. Only a few of the components of 
the Riemann tensor are non-vanishing. These are shown 
at the end of this appendix in subsection A 1. With the 
knowledge of the exact expressions of the Riemann com¬ 
ponents, it is seen that F can be written as: 


2F = TZtrtr (2ie*y{e*rk^k^ - [(e*)l" [k^f - [(e*)f (F)^) + 

Tltptp (2(e*)‘fc‘ {ie*rkP + (e*)«fcn - (fc*)^ {[(e*)^’]' + [(e*)«]'} - [{e*)f {{kP)^ + (A:^)^}) + 

Tlrprp (2{e*y{e*)Pk^kP - [{e*yf {kP f - [(e*)^]" (F)^) + (A15) 

Urqr, (2(e*)^(e*)U«fc’' - [(e*)n" {k'if - [(e*)«]" (F)") + 

Tipqpq {2{e*yy*ykPk‘i - [(e*)^]' (fc«)2 - [(£*)«]' {kpy) 

In order to solve the set of ODEs given above, initial conditions need to be determined. Most of the geodesics are 
initialized as radial. Since = (1,0, 0,0), this implies that the initial conditions for the screen space vectors are 

Ff = ^(0,0,1,0) and El^ = ;^(0,0,0,1). In the more general case of rays that are not initially radial, the screen 
space basis unit vectors are initialized according to: 


Ff (X 



VR klkP^ 

vf ^/{kpy + {kfy 


VR fcgfcg ^ 

vfv/w + w; 


(AI 6 ) 


E^(x 



1 feg 


1 K \ 

Vf ^/Ikpy + V^ J 


(A17) 


These expressions are obtained by requiring that the 
two basis vectors are orthogonal to the velocity field and 
kP. 

In order for y/\ det(F)| to be the angular diameter 


distance, E must be initialized as k^ = — i. The initial 
values of fc'’, fcP, A;'^ are arbitrary though k^ must fulfill 
the null-condition. 

By differentiating the equation ^ it can be seen 

that D must be initialized as the zero-matrix and D as 
the identity matrix. 
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1. Riemann components of Szekeres models 

The Riemann components are computed using the definition TZaPfn, = | [gj)ii,ua — 9afi,i^p + gav,fj.p — gpv,fia) — 
g'^'^ — Tcrau^-fp^i)- This equation contains the covariant Christoffel symbols symmetric in their two last 

indices, i.e.\ 


^ct/3 ~ r,9~^ {ga<7,P + gpa,a 5a/3,cr) 

2 ^ (A18) 

= > r^Q,^ = — {gaj^p + gp~f,a gaP,'Y) 

The non-vanishing covariant Christoffel symbols of the quasi-spherical Szekeres models with the metric corresponding 
to the line element of equation (Al) are: 




r rrt ~ 

■ ^ rtr — 

r trr — 

Rt 

2 





Trrr = 

R,r 



r — 

X - 

P — 

X j.pj, — 

_P — 

X p^y’ — 

to 



r — 

P — 

_P — 

R,q 

2 

II 

II 

1 

II 

= ^qqt = 

= ^qtq = 

-^tqq = 

F,t 

2 

r — r — - 

X ppr ^ prp 

_r — 

X rpp 

P — 

^ qqr 

P — 

^ qrq 

—Y — 

rqq 

F,r 

2 


r — 

1 qqq — 

ppq — 

pqp — 

_P — 

qpp — 

Fq 

2 


r — 

^ ppp — 

P — 

-L qqp — 

P — 

i- qpq — 

_P — 

i- pqq — 

F,p 

2 


The corresponding non-vanishing Riemann component are: 


(A19) 


_ 1 

R-trtr — R-rtrt — R-trrt — Rrttr — '^R,tt T 

_ 1 

R-tptp — Rptpt — R-tppt — Rpttp — T 

_ 1 F2 

R-tqtq — R-qtqt — R-tqqt — R-qttq — 

TZrprp = Rprpr = —Rprrp = —R-rppr = F',rr) + {R,rF,r + R^p) + {R,pF^p + — R,qF^q^ 

R-rqrq = Rqrqr = —Rqrrq = —R-rqqr = {F,qq + T’.rr) + ~r~^F^tF^t + {R,rF,r + R^g) + —^ i^ — R^pF^p + R qFg + 


R'pqpq — R'qpqp — R-pqqp — 'Rqppq 


4A 
1 

4F 

4c2"’' 4R' 


= -5 iF„ + F„) + Tfj _ ^ (fj + fj) 


(A20) 


Appendix B: The Sachs formalism combined with 
the Born approximation 

The ray tracing scheme described in section IV B 
reproduces the exact results very precisely in the special 
case where the exact rays are radial rays in an LTB 
model. Unfortunately, that ray tracing scheme is, in 
practice, difficult to use for studying rays initialized as 


non-radial in the local single void models of the swiss 
cheese. Thus, in order to test if the good reproduction is 
purely a consequence of the exact rays themselves being 
radial, the Sachs formalism has been used to devise a 
simple scheme for ”ray tracing” non-radial rays in LTB 
models. In this ray tracing scheme, rays are propagated 
according to the background null-geodesic equations, 
i.e. in accordance with the Born approximation. 
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The convergence can be obtained using the 
Sachs formalism described in section IV A. Using 


TZ = 




8-irG 


-T 

II 


47rG 


R = = - 


it is easily seen that 
2 Y 2 -p {—cPk* + a^Vrk^)'^. Since 


the ray tracing will only be done through single void 
models, the gravitational convergence is, however, 
negligible. This implies that the result will not depend 
on whether it is the background density or mapped 
Szekeres density that is used when computing R. In 
fact, the term a^Vrk'' can also be neglected, and F is set 
equal to zero, corresponding to its lowest order value. 
The only place where the computations deviate from 
background computations is thus when computing the 
observable redshift which corresponds to computing the 
Doppler converge nce. _ 

In order for yj\ det(II)| to be the angular diameter 
distance, fc* must be initialized as A:q = — ^ while the 
initial values of are arbitrary though the null- 

condition for must be fulfilled. D must be initialized 
as the zero-matrix and D as the identity matrix. 

The formalism described here has only been used 
with fc® initially set to zero. In this case, the ini¬ 
tial conditions of the screen space unit vectors are 
cx (^,^,-A2|^,0) and E!^ = (0,0,0, J^), 
where if is a metric function defined in appendix C. 
These initial conditions are obtained by using the fact 
that the screen space basis vectors are orthogonal to 
each other and to « ^(—c^(l-|-2?/;), a^(l —2^/^)ur, 0, 0) 
and k^. Once again, it is actually sufficient only to keep 
the background terms. 


At low redshifts, the scheme just described has the 
same precision as the ray tracing scheme described 
in section IV B; it uses the Born approximation and 
computes the Doppler convergence. Thus, if the ray 
tracing scheme described here can accurately reproduce 
the path and corresponding convergence of a non-radial 
light ray in an LTB model, then the ray tracing scheme 
of section IV B will also be able to. As shown in section 
V A 1, this is in fact the case. 


rized below together with examples. 


1. Mock N-body data 


Disregarding notation, the map in equation (6) does 
not require the background to be flat. With a gen¬ 
eral background, the right-hand side should refer to the 
given FLRW model of choice and not speciflcally the EdS 
model, i.e.\ 

t .— tsz — ^flrw 
P Psz — Pfirw 

Q •— Qsz — Qflrw (^1) 

rrsz 

I d'^ sz rr,sz — / drfi rw\/f/rr,flrw 

Jo Jo 

The coordinate system used for the FLRW model should 
be the one corresponding to the FLRW limit of the Szek¬ 
eres spacetime. The coordinate choice should thus corre¬ 
spond to the FLRW line element written as: 

+ a> + V)) 

(C2) 

K G {—1,0,1} gives the sign of the spatial curvature 
while i?n = - , is the curvature radius. The 

function E is given by if = ^ -I- -I- l). 

The velocity field cannot be obtained using equa¬ 
tion (7); the right-hand side of that equation is only 
correct for a flat FLRW model. The velocity held can be 
obtained through the same considerations as those that 
lead to equation (7): 

The point identification map of equation (Cl) is time 
dependent. Hence, keeping the spatial Szekeres position 
constant while letting the time coordinate vary will lead 
to a mapping into spatially different FLRW points. This 
change in the spatial FLRW coordinates as a function of 
the time coordinate leads to a velocity field that can be 
computed e.g. by using a finite difference formula. 


Appendix C: Curved FLRW backgrounds 


2. Ray tracing 


The truncated onion models shown in figure 2 all con¬ 
verge to an EdS background at large r. The structures 
at smaller values of the radial coordinate are however 
not separated from each other by artificial background 
patches but are free to interact. This implies that the 
structures will generally not be locally compensating. 
Thus, ray tracing through the inhomogeneous region of 
the onion models may require the use of a different FLRW 
model than the EdS model as the background. Some of 
the equations of sections III and IV B require the back¬ 
ground FLRW model to be flat and thus need to be 
adapted before they can be used together with a curved 
background. The necessary modifications are summa- 


The expression for the gravitational convergence for a 
general FLRW background is well known (see e.g. [74, 
84]): 


KS = 


47rG 


dx'a^Sp 


2. fK{Xg-x'g)fK{x'g) 


fKiXg) 


(C3) 


The relation between x and the radial coordinate is given 
by r = fK{x), where fxix) is the curvature function 
given by: 


I i?o sin(x/i?o) if AT = -1-1 
fK{x) = \ xiiK = 0 (C4) 

I Rq sinh(a;/i?o) if AT = —1 
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Density along radial ray, background with = 0.8 convergence along radial ray, background with Qrn,o — 0.8 



^obs 


0.05 


(a) 



FIG. 16. Density and convergence along a radial geodesic in an LTB swiss cheese model. The density profiles along the exact 
and ray traced rays are indistinguishable. The same is true for the exact convergence and the sum of the gravitational and 
Doppler convergence. 


Equation (C3) requires the scale factor to be normalized 
at the time of observation in the FLRW background co¬ 
ordinates. According to the map in equation (Cl), this 
time coordinate is the same as the time coordinate of 
the observer in the exact Szekeres spacetime. The Hub¬ 
ble constant of a given background must thus be chosen 
such that the age of the background model is the same 
as that of the exact Szekeres model being studied. 

The expression for the Doppler convergence computed 
on a curved FLRW background can be written as: 


KyC = i— 


^n • (v - Vo) -f n • Vo 

(C5) 


The radial coordinate r is the coordinate corresponding 
to a point along a radial null-geodesic in the FLRW 
background at z = Zobs- 

The expression can e.g. be obtained by following 
the procedure presented in section 2.2 of [99], modified 
to a curved background. The main consequences of 
introducing a curved background are due to changes 
in the expression for the background angular diameter 
distance and in particular its dependence on the redshift. 

For the studies presented here, the factor 
- Kr^/Rl is approximately equal to 1 and may 
thus be neglected both in the expression for the Doppler 
convergence and in the relation between r and the 
conformal time ry. Consequently, the expression for the 
Doppler convergence shown above may be approximated 
by the expression for the Doppler convergence on a flat 
FLRW background. 


The ray tracing method on a curved background 
has been tested using different swiss cheese models. An 
example is shown in figure 16. The figure shows the 
density and convergence along a ray moving radially 
through consecutive LTB voids. The particular single 


void model used to construct the swiss cheese model has 
been specified similarly to the models in section II but 
with k{r) given by: 



The results depicted in figure 16 correspond to setting 
= 0.8. This gives a model that converges asymp¬ 
totically to an open FLRW model with a present day 
density parameter of = 0.8. The ray tracing results 
have been obtained by using this as the background 
model. As indicated in the figure, the ray tracing results 
reproduce the exact results very precisely. 

The model defined by equation (C6) has an unrealistic, 
pointy void shape indicating that the model is not well 
behaved at the origin. Since the results are reproduced 
well by the ray tracing scheme this is not an issue for 
the current purpose of the model. The large patches of 
nearly homogeneous background density along the ray 
are necessary to avoid velocity artifacts. 

When changing the background in the ray tracing 
scheme, the background angular diameter distance 
computed following equation (10) and used in equation 
(9) must also be changed accordingly. 


3. Ray tracing through swiss cheese models with 
improper backgrounds 


As shown in section V A, the ray tracing scheme very 
precisely reproduces the exact convergence along rays 
moving radially through consecutive LTB voids. It is 
interesting to see to what extent this result requires using 
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the proper background. By ’’proper background” is here 
meant the FLRW patches between the compensated 
single void structures of the given swiss cheese model. 
For the study in section V A, the proper background 
was the EdS model. Figure 17 shows the convergence 
obtained by using different curved backgrounds for the 
ray tracing approximation and in equation (9). The 
Swiss cheese model is the same as the one studied in 
section V A. 

As seen in figure 17, the background is important 
in order to obtain precise results with the ray tracing 
scheme. It is as expected that the ray tracing should 
yield an incorrect result when using an improper back¬ 
ground; the Doppler convergence becomes incorrect 
when computed from an improper background since 
the peculiar velocity field is supposed to describe the 
difference between the local exact expansion of the LTB 
model and the global expansion of the average model. 
This is only the case if it is the proper background that 
is used. Note in this connection that when improper 
backgrounds are used, the mock N-body velocity field 
does not reduce to zero outside the single void structures 
of the cheese. 

Using an improper background for computing the 
exact convergence from equation (9) should lead to a 
shift so that the exact convergence no longer oscillates 
about zero. This is clearly seen in figure 17. This 
same does not happen in the ray tracing scheme where 
the Doppler convergence keeps oscillating about zero. 
While the gravitational convergence is negligible, the 
ray tracing convergence will thus oscillate about zero 
regardless of the background. 

The LTB onion models have also been studied us¬ 
ing curved backgrounds. This was done as an attempt to 
obtain a better reproduction of the exact convergence. 
This turned out not to be possible. In fact, the differ¬ 
ence between exact and ray traced convergence does not 
appear to depend much on the used background; the 
exact and approximate convergences are simply raised 
or lowered by approximately the same amount when 
the background is changed. This is very different from 
the Swiss cheese case illustrated in figure 17 where the 
difference between the exact and ray tracing convergence 
clearly depends on the background. The cause for 
this could be related to the very different Doppler 
convergence of the onion models compared to that of 
the LTB Swiss cheese models. 
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convergence along radial ray, background with ^^,0 = 0.7 



{^)^m,o — 0-^ 


convergence along radial ray, background with 0,^,0 = 0.95 



{c)nm ,0 = 0.95 


convergence along radial ray, background with = 1-01 



(e)Slm ,0 = 1-01 


convergence along radial ray, background with q = 1-1 



(s)^TTT',0 — 1-1- 


convergence along radial ray, background with = 0.9 



(b)Qm,o = 0.9 


convergence along radial ray, background with = 0.99 



(d)nm,Q = 0.99 


convergence along radial ray, background with = 1-05 



(f)!^m,o = 1.05 


convergence along radial ray, background with Qrn,Q = 1-3 



(h)Q^,o = 1.3 


FIG. 17. Convergence along a radial ray in the LTB swiss cheese model based on the LTB single void model specified by m = 4 
and Vb = 60Mpc. The exact and ray traced convergences are computed using various FLRW backgrounds specified by their 
values of flm.o- The proper background is the EdS model. 
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